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Abstract. We study the three nucleon force in Nt = 2 dynamical clover fermion lattice QCD, 
utilizing the Nambu-Bethe-Salpeter wave function of the three nucleon system. Since parity-odd two 
nucleon potentials are not available in lattice QCD at this moment, we develop a new formulation 
to extract the genuine three nucleon force which requires only the information of parity-even 
two nucleon potentials. In order to handle the extremely expensive calculation cost, we consider 
a specific three-dimensional coordinate configuration for the three nucleons. We find that the 
linear setup is advantageous, where nucleons are aligned linearly with equal spacings. The lattice 
calculation is performed with 16 3 x 32 configurations at j3 = 1.95, m n = 1.13GeV generated by 
CP-PACS Collaboration, and the result of the three nucleon force in triton channel is presented. 
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INTRODUCTION 

Understanding the interaction among nucleons is one of the most important issues in 
, particle and nuclear physics. In particular, the properties of the three nucleon force 

(TNF) are the forefront of not only particle/nuclear physics but also astro physics. For 

^ ■ the binding energies of light nuclei, the two nucleon (2N) force alone underestimates the 

experimental data, and attractive TNF components are required [1]. On the other hand, 

1/-^ ; repulsive TNF components are also necessary in both of symmetric and asymmetric 

nuclear matters, in order to reproduce the empirical saturation density of nuclei in the 
former case, and to make the equation of state of asymmetric nuclear matter hard enough 
to sustain the observed massive neutron star masses in the latter case [2]. Recently, it is 

. £h 1 argued that TNF is responsible for the anomaly in the drip line and the nontrivial magic 

number of neutron-rich nuclei [3, 4], and could impact on the supernova nucleosynthesis. 

Pioneered by Fujita-Miyazawa [5], TNF have been commonly studied from the two- 
pion exchange models with the A-excitation. However, since TNF is originated by the 
fact that a nucleon is not a fundamental particle, it is desirable to study TNF from the 
fundamental degrees of freedom (DoF), i.e., quarks and gluons. In this report, we carry 
out first-principle calculations of TNF using the lattice QCD technique. 

In the lattice QCD calculation of nuclear forces, it is recently proposed [6] to use 
the Nambu-Bethe-Salpeter (NBS) wave function so that the potential is faithful to the 
phase shift by construction. The obtained (parity-even) 2N potentials are found to have 
desirable features, such as attractive well at long and medium distances, and the central 
repulsive core at short distance [6, 7]. It is also examined that the non-locality of the 
obtained potential is small, and the leading terms in the derivative expansion dominate 
the potential at low energies [8]. The method has been successfully extended to the 
hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions [9, 10, 11] including 



H-dibaryon system [12]. Further applications are given in Refs. [13, 14, 15, 16]. In 
this report, we discuss the extension of this approach to the three nucleon (3N) system, 
and lattice QCD calculation of TNF in the triton channel is carried out. A preliminary 
account of this study is also given in Ref. [17]. 



FORMALISM 

Since the detailed formulation for the 2N force is given in Ref. [7], we discuss the 
extension to the 3N system. We consider the equal-time NBS wave function y(r , p), 
which can be obtained through the six point correlator, 



G a py,aipif{r',P,t-to) = (N a (xi,t)N p (x 2 ,t)N r (x3,t) (A^(? )%(A))AV (?()))>, (1) 

where r' = x \ — X2, p = xj, — (x\ +X2)/2 are the Jacobi coordinates, and N denotes either 
of p or n. With the derivative expansion of the potentials, the NBS wave function can be 
converted to the potentials through the following Schrodinger equation, 
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\lf(r J ,p)=E\jf(r',p), (2) 



where V2N,ij{xij) withjc/y =Xi—xj denotes the potential between (/,y')-pair, Vtnf^ \P) 
the TNF, /v = «%/2, p p = 2m^/3 the reduced masses. If we calculate y(r , p) for all ? 
and p, and if all V2N,ij(xij) are available by (separate) lattice calculations for the genuine 
2N system, we can extract VrNFij 1 ,p) through Eq. (2). 

In practice, however, the calculation cost of \\f{r' ',p) is prohibitively expensive. Note 
that the number of diagrams generated by Wick contractions tends to diverge with a 
factor of N u l x Njl (N Uj d are numbers of u,d quarks in the system). We develop several 
techniques to reduce the calculation cost, e.g., we take advantage of symmetries (such as 
isospin symmetry), and employ the non-relativistic limit for the source nucleon operator. 
Similar techniques are (independently) developed in the calculations of binding energies 
of triton/helium nuclei [18]. Yet, the calculation cost remains quite expensive, so we 
consider Eq. (2) with fixed three-dimensional (3D) coordinates of r' and p. 

To proceed, we find that the linear setup with p = is advantageous, where three 
nucleons are aligned linearly with equal spacings of r = \r\ = l/2|r'|. In fact, because of 
p = 0, the third nucleon is attached to ( 1 , 2) -nucleon pair with only S-wave. Considering 
the total 3N quantum numbers of / = 1/2, J p = 1/2+ (triton channel), the wave function 
can be completely spanned by only three bases, which can be labeled by the quantum 
numbers of (l,2)-pair as ^o, 3 ^i, 3 £>i. Therefore, the Schrodinger equation leads to the 
3x3 coupled channel equations with the bases of s , y/3 5 , . The reduction of the 
dimension of bases is expected to improve the S/N as well. 

Even when we calculate with the fixed 3D-configuration (with the linear setup or 
whatever), the identification of genuine TNF remains nontrivial. This is because the 
parity-odd 2N potentials are not available in lattice QCD at this moment, and we cannot 
subtract the contribution of V^s unambiguously. Note that although the total parity of the 



3N system can be projected, a 2N-pair inside the 3N system could be either of positive or 
negative parity. The familiar procedure of partial wave expansion cannot be performed 
here, since we calculate only fixed 3D-configurations to reduce the calculation cost. 
To overcome this problem, we find that the following channel in the triton is useful, 



This wave function itself has been well known, but the point here is that it is anti- 
symmetric in spin/isospin spaces for any 2N-pair. Combined with the Pauli-principle, 
it is automatically guaranteed that any 2N-pair couples with even parity only. Therefore, 
we can extract TNF unambiguously in this channel, without the information of parity- 
odd 2N potentials. Note that no assumption on the choice of 3D-configuration is imposed 
in this argument, and we can take advantage of this feature in future TNF calculations 
with 3D-configurations other than the linear setup. 



LATTICE QCD RESULTS 

We employ Nf = 2 dynamical clover fermion configurations generated by CP-PACS 
Collaboration [19]. We use 598 configurations at j8 = 1.95 with the lattice size of 
L 3 x T = 16 3 x 32, which corresponds to (2.5fm) 3 box in physical spacial size with 
the lattice spacing of cT x = 1.269GeV. We calculate at the hopping parameter of u, d 
quarks k uc i = 0.13750, which corresponds to m n = 1.13 GeV, = 2.15 GeV. We use 
the wall quark source with Coulomb gauge fixing. In order to enhance the statistics, we 
perform the calculation for sources on 16 time slices for each configuration. 

We perform the calculation with the linear setup at ten physical points of the distance 
r. Since we have only one channel (Eq. (3)) which is free from parity-odd 2N potentials, 
we can determine one type of TNF. In this report, we consider the scalar/isoscalar 
type TNF. In the Urbana IX model [20], which is an often-used phenomenological 
TNF, TNF consists of two parts, one is the two-pion exchange TNF and the other 
is phenomenologically introduced scalar/isoscalar repulsive TNF. Since the two-pion 
exchange TNF is expected to be suppressed in the current lattice setup of m % = 1.13 
GeV, it is reasonable to assume the scalar/isoscalar nature for TNF here. 

In Fig. 1 (left), we plot each wave function of \j/s = (— Vis + V / 35 1 )/v / 2, Wm = 
(V^Sq + V^sJ/v^ and ty3 Dl in the triton channel at t — to = 8. We observe that \ys 
dominates the wave function. This is because Xj/g contains the component for which 
all three nucleons are in S-wave. By solving the 3 x 3 coupled channels and subtracting 
the V2N from the total potentials, we determine TNF. In Fig. 1 (right), we plot the results 
for the scalar/isoscalar TNF. Here, the r-independent shift by energies is not included, 
and thus about £?(10) MeV systematic error is understood. At the long distance region 
of r, TNF is small as is expected. At the short distance region, we observe the indication 
of repulsive TNF. Recalling that the repulsive short-range TNF is phenomenologically 
required in high density matter, this is very encouraging result. Of course, we note that 
further researches are necessary, e.g., the examination of the discretization error and the 
dependence on the quark masses. Studies along this line is currently underway. 
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FIGURE 1. (color online). (Left) The wave function with linear setup in the triton channel. Red, blue, 
brown points correspond to ty/j, \]/m, V^d,' respectively. (Right) The scalar/isoscalar TNF in the triton 
channel, plotted against the distance r in the linear setup. 
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